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Abstract 

Multiuser Multi-Packet Transmission (MPT) from an Access Point (AP) equipped with multiple antennas to 
multiple single antenna nodes can be achieved by exploiting the spatial dimension of the channel. In this paper 
we present a queueing model to analytically study such systems from the link-layer perspective, in presence of 
random packet arrivals and heterogeneous channel conditions. The analysis relies on a blind estimation of the 
number of different destinations among the packets waiting in the queue, which allows for building a simple, but 
general model for MPT systems. Simulation results validate the accuracy of the analytical model and provide further 
insights on the cross-relations between the channel state, the number of antennas, and the number of active users, 
as well as how they affect the system performance. The simplicity and accuracy of the model makes it suitable for 
the evaluation of link-layer protocols supporting multiuser MPT in non-saturation conditions, where the queueing 
dynamics play an important role on the achieved performance. 

Index Terms 

MPT, SDMA, MIMO, Link-layer, Queueing model, Wireless Networks 

I. Introduction 

In packet-based wireless networks, the use of spatial multiplexing allows for simultaneous transmission 



of multiple packets, directed to a single or multiple destinations, hence increasing the system capacity [1]. 
_C3 We focus on a scenario where an AP equipped with multiple antennas is able to simultaneously transmit 
multiple packets, each directed to a different single-antenna node. This scenario is known as Multiuser 
Multi-Packet Transmission (MPT), a packet-based extension of the Multiuser MIMO downlink [2], also 
referred to as Space Division Multiple Access (SDMA), depending on the context. 

Research on Multiuser MPT has mainly been focused on the design of efficient joint precoding and 
user selection strategies, as a trade-off between computational complexity and the ability to maximize 
the system capacity, i.e., the number of bits per Hertz of available bandwidth that can be successfully 
transmitted over the channel. A comprehensive survey of such works has been presented in [3]. In these 
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works, it is usually assumed that the transmitter has separate per-user queues that are always saturated. 
This reduces the problem to finding the set of users that maximize the system sum-rate capacity, based 
on the current state of the channel. Specific schemes for user selection range from random selection to 
greedy schedulers that benefit from the existing multiuser diversity [4]. 

In contrast, there are still very few works that address spatial multiplexing from the link-layer 
perspective, and even fewer that focus on the queueing dynamics. In [5], the authors present a point- 
to-point MIMO system, considering both transmit diversity (STBC) and spatial multiplexing (BLAST) 
schemes for a single-user MPT. A similar work that overcomes the approximation done in [5] for the 
BLAST scenario is presented in [6], although the presented results are only valid for two antennas. 

To the best of our knowledge, [7] is the only work where a queuing model for modeling the Multiuser 
MPT is presented. In their work, the nodes are assumed to be equipped with at least as many antennas 
as the AP. Therefore, the AP can use its multiple antennas to send data to a single user or multiple users 
simultaneously, following the independent stream scheduler approach [8]. Specifically, they consider the 
V-BLAST scheme [9], using a Zero Forcing (ZF) detector at the nodes, which are also equipped with 
multiple antennas. 

In this work we present a queueing model for Multiuser MPT systems following a completely different 
approach compared with [7], where the authors focus on the analysis of a single tagged user, by setting 
our focus on the complete system. Moreover, in our system model, the nodes are assumed to have only 
a single antenna, and therefore, the AP cannot transmit more than one packet to each node at any given 
transmission. Therefore, the number of packets included in each transmission depends on the number 
of nodes represented in the queue at the moment the transmission is scheduled. Our analytical model is 
based on the blind estimation of the number of different destinations among the packets waiting in the 
queue [10], which allows us to build a computationally simple but general and accurate queueing model 
for Multiuser MPT systems. Therefore, the presented model can be used to understand and evaluate the 
different interactions that exist in a Multiuser MPT system between the traffic load, the buffer size, the 
number of antennas at the AP, the number of different nodes, the channel characteristics, and the protocol 
overheads required for CSI estimation and reporting. In addition, due to the model characteristics, it can 
be easily coupled with other link-layer mechanisms to evaluate more complex systems. 

The paper is structured as follows. The scenario, together with the system model and the assumptions 
considered, is introduced in Section II. Section III presents the Multiuser MPT queuing model. Section 
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IV presents the results, including the validation of the queueing model. Finally, the main conclusions of 
this work are summarized and the future research lines are stated. 



II. System Model and Assumptions 



A network consisting of a multiple- antenna AP and N single- antenna nodes located a single hop away 
from the AP is considered. The AP is equipped with M antennas and M radio-frequency (RF) chains, 
which allow it to transmit up to M simultaneous packets. Packets of a constant length of Ld bits destined 
to the N single-antenna nodes arrive to the AP according to an aggregate Poisson arrival process of rate 
A, containing independent and identically distributed shares of traffic per node. The AP has a finite-buffer 
of size K packets. All nodes have the same priority, and therefore, the finite-buffer space is fully shared 
by all arriving packets, which is optimal if no traffic differentiation is required [11]. A detailed model of 
the AP architecture is shown in Figure 1. 

A. Space-Batch Scheduling 

At every transmission opportunity, the AP starts by constructing what we refer to as a space-batch, 
containing up to M packets. The space-batch is constructed on a First-Come-First-Serve (FCFS) basis, 
however, once a packet destined to a certain node is placed in the space-batch, subsequent packets to the 
same destination are skipped and left for future space-batches. This is because each node is only equipped 
with a single antenna and cannot receive more than one packet at a time, and therefore, there can be at 
most one packet per destination in every space-batch. 

The scheduling of space-batches takes place immediately after the completion of each transmission if 
the queue is not empty. Otherwise, the AP will wait until a new arrival enters the queue, immediately after 
which a space-batch containing only one packet will be constructed. When there are strictly more than 
one packets present in the queue after a departure, there may be multiple packets among them destined 
for the same node. Therefore, only a subset of the packets in the queue might be eligible for transmission 
in a single space batch. Let £ be the number of eligible packets (destined to distinct destinations) in the 
queue immediately after a given departure. Then the size of the space-batch scheduled after that departure 
is given by: 

1 £ < 1 

s(£) = max(l, min(£, s max )) = <J £ 1 < £ < s max (1) 

$max Smax ^ £ 
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where s rnax < M is a system constant indicating the maximum number of spatial streams allowed to be 
sent at each transmission. Choosing a value for s max that is strictly smaller than M enables the transmitter 
to benefit from extra transmit diversity gain, which in turn results in higher SNR at the receiver. 

B. Frame Structure 

Let S be the set of nodes for which a packet is included in the current space-batch. After constructing 
the space-batch, the AP broadcasts the identity of the nodes in S. This is immediately followed by a 
series of training sequences from the AP to the selected destinations. Based on the training sequences, 
each node in S reports its estimated CSI back to the AP. The AP will use the received CSI to form the 
required beams and to choose the appropriate transmission rate, as will be detailed shortly. 

The general structure of a transmitted frame is depicted in Figure 2. It consist of five parts, transmitted 
in the order presented below: 

• Preamble and space-batch information (of length L s & bits): Contains the initial preamble used to 
synchronize all receivers and the headers required to inform those nodes that have been selected to 
receive a packet in next space-batch. 

• Training Sequences (of length M x L tr bits): Required to estimate the channel state between each 
transmitting antenna at the AP and the receiving antenna at each selected node. 

• CSI feed-back (of length m x L C si bits): Used for nodes to report their estimated CSI to the AP. 

• Data Packet (of length L d bits): Includes the space-batch data packets. 

• ACKs (of length m x Lack bits): Used for nodes to notify the correct reception of the data packet. 

C. Channel Model and Transmission Rate Selection 

A quasi-static fading channel that changes from space-batch to space-batch transmission is considered. 
Let r n m denote the instantaneous Signal-to-Noise ratio (SNR) observed by node n E S when m spacial 
streams are included in the space-batch, i.e., m = \S\. The instantaneous SNR observed by a node is 
assumed to be independent of that of the other nodes. Considering that a ZF beamforming is used, the 
received SNR at each node is assumed to follow a ^-distribution with / = 2 x (M — m + 1) degrees of 
freedom [12]. The cumulative distribution function for r nm is therefore given by: 
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where T n . m is the average SNR observed by node n when a space-batch of size m is transmitted and 
depends only on the pathloss between the AP and node n. Note that f n>m = T n> i/m, since the transmission 
power is equally divided between the m parallel streams. 

Upon receiving the CSI feedback, which is assumed to be ideal, the AP picks a transmission rate based 
on the SNR values calculated for every selected node. The transmission rates are chosen from a finite set 
of values 1Z = {r ,ri, . . . ,tr}, where r = < r\ < ■ ■ ■ < Tr. When r is selected, no transmission is 
possible. The transmission rate for the n-th. spatial stream, f n , is chosen to be if T n m falls in the range 
(7i,7i+i], where {^j}f=Q 3X6 predetermined thresholds. Once the transmission rate at which each spatial 
stream can be transmitted is known, the AP chooses the smallest rate, f = mm{f n } n€S , and uses it as 
the transmission rate for the whole space-batch. 

It is assumed that the probability of a space-batch suffering transmission errors is negligible if the 
proper f is used. In the specific case where f = r , the channel state does not allow for the transmission 
of the current space-batch, and the AP has to defer any transmission attempt for a period equal to T a , in 
hope of a better channel state. After this time period, a new space-batch will be scheduled based on the 
current state of the queue and a new frame will start accordingly. From here on, T Q is referred to as the 
channel outage period. 

It should be noted here that the control parts of the frame are all transmitted omni-directionally at r\, 
the smallest non-zero rate, in order to maximize the probability that all nodes in the coverage area of the 
AP can successfully decode them. The duration of the control part of the frame is independent of the 
chosen transmission rate, r, and is given by: 



The duration of a frame depends on the number of antennas and the transmission rate r chosen for the 
space-batch transmission, hence: 



D. An Example of the System Operation 

In Figure 3, a specific example of the system operation is shown for M = 2 antennas and a queue 
of K = 4 packets in a system with transmission rates 1Z = {r , r±, r 2 , r^}. The (i — l)-st space-batch 



T c (m) 



L sb + M ■ L tr + m ■ L CS i + m ■ L ACK 



(2) 




(3) 
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comprises a single packet as the transmission is scheduled as soon as a new packet arrives to the AP and 
it can be transmitted at rate r 3 . During the (i — l)-st space-batch transmission, two packets, one directed 
to the 4-th node and the other to the 3-rd node, arrive to the AP and are buffered and assembled together 
in the i-th space-batch after the (i — l)-st space-batch transmission is completed. However, in this case, 
the channel does not allow the transmission of the scheduled space-batch (f = r ) and therefore, the AP 
decides to defer the transmission and schedule a new space-batch after T a seconds. During this period 
another packet directed to the 4-th node arrives to the AP and is placed in the queue, as it still has free 
space available. After the T a seconds, a new space-batch is scheduled containing a packet each to nodes 
3 and 4, which at this time can support rates f 3 = r 3 and r± — r\, respectively. The space-batch is then 
transmitted at rate r — r\ — min{r3,ri}. Similarly, during the z-th space-batch transmission two more 
packets directed to the 4-th node arrive to the queue, as well as one directed to the 5-th node, which 
is blocked because there is no free space in the buffer. Observe that, when the [i + l)-st space-batch is 
scheduled, there are only two packets in the transmission buffer and both are directed to the 4-th node. 
In such situation, only one packet will be transmitted, which in our example is done at rate r 2 . 

E. Discussion of System and Modeling Assumptions 

In the previous described system model, there are several assumptions which require further justification 
as they are done to make possible the model, or have implications on its performance. Nevertheless, they 
are general enough to provide a clear understanding of the system operation, which can be used as a 
reference for future developments. 

• Poisson Arrivals and Homogeneous traffic: The assumption of Poisson packet arrivals has been 
traditionally used in queueing models due to its specific properties [13]. One of those properties, 
known as Poisson Arrivals See Time Averages (PASTA) property, is used in this paper to simplify 
the computation of the probability distribution at arbitrary times as, for the case of Poisson arrivals, it 
is equivalent to the distribution at packet arrivals. Regarding the assumption of homogeneous traffic, 
it allows us to find a closed-expression for the blind estimation of the space-batch size distribution, 
which otherwise may not be possible. 

• Blind estimation of the space-batch size: As we will see in Section III, the model is build on top 
of the assumption that we can estimate the distribution of the space-batch size only based on traffic 
distribution among nodes and the number of packets waiting for transmission when a new space-batch 
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is scheduled. In the one hand, this approximation allows us to build a simple but accurate model, as 
it will be shown in Section IV although, in the other hand, it has an implicit limitation. It is that the 
model does not known the specific destinations of the selected packets, which prevents any further 
action once a certain number of packets have been selected for transmission. 

• Outage periods: The fact that m packets could not be transmitted together (i.e. f = r ) does not 
mean that less than m packets cannot be successfully transmitted. However, to solve this point, the 
AP has to search if there is any suitable combination between the m selected packets, which means 
that it has to know what are the selected destinations. As discussed in previous point, this is not 
straightforward following our modelling strategy. 

• Transmitting the space-batch at the smallest r: As it is assumed that all packets have the same 
length, the duration of a space-batch is always limited by the smallest transmission rate. In such 
case, transmitting the other packets at higher rates does not result in any performance gain or change 
in the system dynamics as no more packets directed to the nodes with higher rates can be scheduled 
until the overall space-batch transmission finishes. To improve the system performance, a solution 
can be to consider packet aggregation / fragmentation in order that given a certain duration, each 
node could transmit a different amount of bits depending on its specific rate. 

III. QUEUEING IN MPT SYSTEMS 

A. Model 

An M/G^/l/K queue is used to model the AP. In this notation, [s] = [1, s max ] indicates an interval, 
with 1 and s max representing the minimum and maximum number of spatial streams that can be used 
at each transmission, respectively. Additionally, Poisson arrivals of rate A and a general service time 
distribution are considered. The buffer has a size of K packets and no extra space is considered for the 
packets in service, i.e., the packets included in a space-batch transmission remain stored in the queue 
until their transmission is completed. 

As explained before, the scheduling of a space-batch takes place immediately after each departure, if 
the queue is not empty. Here a departure instant refers both to the moments at which all packets from a 
successfully transmitted space-batch are purged from the queue and to the end of channel outage periods. 
If just after a departure the queue is empty, the AP will wait until there is a new arrival, at which instant 
a new transmission containing a single packet is scheduled. 
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Let q k denote the queue occupancy at the end of the A;-th interdeparture epoch. If m k is the number 
of successfully transmitted packets and r(k) the transmission rate at which the space-batch is transmitted 
during the A;-th epoch, then q k evolves according to the following recursion: 

q k = min {v[T(m k , f{k))\ + -m k ,K - m k } (4) 

where v[T(m k ,r(k))} is the number of packet arrivals during the k-th inter-departure epoch. Therefore, 
v[T(m k ,f(k))] + q k _\ — m k is the number of packets that would be present in the queue at the end of 
the k-th. epoch if the buffer had infinity capacity, and K — m k is the maximum possible queue occupancy 
immediately after a departure for a finite buffer. Note that in this notation, m k = when the selected 
transmission rate is f(k) = r . 

In order to find analytical formulation for key performance metrics such as delay and throughput, the 
steady-state queue occupancy probabilities need to be calculated. To derive these probabilities, first, the 
steady-state distribution for the queue states immediately after departures, 7r d , is derived using a discrete- 
time embedded Markov chain. Then the PASTA (Poisson Arrivals See Time Averages) property of the 
Poisson arrivals is applied to find the occupancy distribution at arbitrary times, ir s , as a function of n d . 

In what follows, in order to make a clear distinction between the aforementioned two different steady- 
state probabilities, we define two different sets of states and a corresponding terminology for their 
probability distribution as follows: 



. queue state at departure instants: number of packets stored in the queue immediately after a 
departure. Hereafter, the steady-state probability distribution for these states is referred to as the 
departure distribution. This is what was denoted above by the row vector 7r d . 

• queue state at arbitrary times: number of packets stored in the queue at any arbitrary time. The 
steady-state probability distribution for these states will simply be referred to as the steady-state 
distribution of the queue. This is what was denoted above by the row vector tt s . 



1 ) Blind Estimation of the Space-Batch Size Distribution: Let p m \ q be the probability that m packets 
are scheduled given that q packets are stored in the queue. This probability is given by 



P 



m\q 



1, m — 1, q = 

Pr{5 q = m}, m < s max , q > m 

J2t= Smax Pr { H ? = £}> m = s max, q>m 

0, otherwise 



(5) 



where E q is the random variable denoting the number of packets eligible for transmission in a single 
space-batch, among the q packets present in the queue. 

Each of the q packets in the queue can be destined to any of the N different nodes in the network, 
and therefore there are a total of N g different possible queue arrangements in terms of destination 
representation. To compute the probability Pr{S g = £}, we need to find the fraction of such arrangements 
for which there are exactly £ different destinations represented. 

Consider one such favorable arrangement with exactly £ nodes represented. Let fi^ = (/xi, . . . , /ig) be 
a vector containing the number of packets of each of the £ destinations represented in the queue. Then 
Hi need to be positive integers and Yl\=i A*« = Q- We define ^ to be the set containing all such vectors, 
i.e., ty q £ = {p,£ G | Yli=i = <?}• Then corresponding to every vector in ^ q ^, we have PR^ e possible 
different queue arrangements, where 

PR? = (6) 

is the number of permutations of q elements, partitioned into sets of /xi, . . . ,/xg repeated elements. The 
probability of having £ eligible packets in the queue is then given by: 

(N\ pR g 

Pr{S, = {}= ^W«« *i (7) 

where X^e^ 9 ? 1S tne tota ^ num ber of possible queue arrangements for a fixed set of £ represented 
nodes and (^) is the number of such sets. 

It should be noted that here that in our calculations above it is assumed that, given a randomly chosen 
packet from the queue, the probability that it is destined to any given target node is \/N, i.e., the packets 
in the queue at departure instants are equally distributed among all nodes, in the same way that the 
aggregate traffic is uniformly distributed among all nodes. However, this is not exactly true because the 
space-batches are constructed containing no repeated packets and every space-batch departure will reduce 
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the diversity of the queue. This assumption greatly simplifies the analysis, however, as we will see, it 
does not bear any significant impact on the analytical results, which actually match the simulations quite 
well. 

Equation (7) generalizes the results in [14] by providing a single expression for any value of M, N 
and q. In [14] only a single case is provided for q = 4, N = 4 and M = 4. In this case, since M = N, 
Pm\q — P r {^q = m } for Q > 0, which for q = 4 and different values of m is given by: 



4 /4!\ 

G) / 4! 4! 4! 



P214 = 4^ (vl!3! + 2!2! + 3!lT J (9) 
(3) / 4! 4! 4! 



Pm ~ 4^ VTIT!2! + T!2!lT + 2!1!T! ^ 
1 / 4! 



2) Distribution of the Selected Rate, f: Based on the rate selection mechanism explained in Section 
II-C and the SNR distribution indicated by (2), the probability that a given node n has a feasible rate 
is given by 

9 nAm = Pr{r n = n\m} = F Tnm (y l+l ) - F Tnm {^i) (12) 

where m is the number of packets included in the space-batch. Given the set of destinations, S, included 
in the next space-batch transmission, a rate f = n will be selected if, based on the CSI feedback from 
those nodes in S, the smallest feasible rate is r^, i.e., min{f n } n6 s = r*j. Therefore, <p ri {S), the probability 
that for a given space-batch composition, S, the rate is chosen, is given by: 

4> U {S) = Pr {min{f n } ne5 = n} 

= J Pr{f n > Ti I m} - JJ Pr{f„ > r i+1 \ m} 

= n(E^)-n(E^>) (i3) 

where the second equality is due to the assumed independence of the observed SNR values for different 
nodes and can be interpreted as the probability that all the nodes in S have rates no smaller than T{ but 
not all of those rates are strictly larger than rv 
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3) Distribution at departure epochs, ir d : The departure probability distribution, n d , is computed by 
solving the linear system: 

ir d = ir d P (14) 

together with the normalization condition: 

7T d l T = 1 (15) 

where P is the probability transition matrix of the embedded discrete-time Markov chain of the occupancy 
of the batch-service queue, immediately after departure instants, with each element pij, i,j e [0, if], 
representing the probability to move from state i to state j. 

In this chain, transitions occur at departure instants, i.e., immediately after the complete transmission of 
a frame or at the end of a channel outage period, and the states represent the queue occupancy immediately 
after a departure. The p^j transition probabilities can be viewed as 

Pij = Pr {q k = j\q k -i = i} (16) 

for any k, where q k and q k ^i are related according to the recursion in (4). As it can be seen in the 
recursion, p it j not only depends on the size of the transmitted space-batch, but also on the time it takes 
to transmit the corresponding frame, which depends on the chosen rate f, which in turn depends on S, 
the composition of the space-batch. 

Let pij(m,r) denote the conditional transition probability given that at state i a space-batch of size m 
packets is transmitted at rate r. This probability depends on the value of V[T(m, r)], the random variable 
representing the number of arrivals during the transmission of the m packets sent in state i at rate r. For 
any state i in the chain, the last reachable state is j = K — m. Therefore, 

{Pr {V[T(m, r)] = j — (i — m)} , j < K — m 
(17) 
Pr {V[T(m, r)] > j — (i — m)} , j = K — m 

where i G [0, if], j > i — m, and m < s(i), with the function s(-) as defined in (1). For all other 
values of we have Pij(m, r) = 0. Note that departing at state j = if — m means that the queue has 
been containing q = if packets just before the departure, and therefore, some arrivals have possibly been 
blocked. For all other reachable states from state i, the queue has had room for more packets just before 
the departure and therefore no arrivals could have been blocked. 
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Note that in case of a channel outage, i.e., when r = r , we need to set that no packets are transmitted 
on the right-hand- side of (17) for it to work correctly, i.e, 



Pr{V[T(m,0)]=3-i}, J<K 
Pi,j{m,r ) = ^ (18) 

Pr{V[T(m,0)]>j-i}, j = K 

where T(m, r ) = T c + T Q , as indicated by (3). 

For Poisson arrivals of rate A, the number of arrivals during T(m, r) has in general the following 
distribution: 

poo ( \t) v 

Pi{V[T(m,r)} = v}= / e ~ xtK -^-f T{mtr) {t)dt (19) 

Jo v - 

which in our case, as given m and r the frame duration T(m, r) is constant, it can be simplified to: 

Pr {V[T(m,r)} = v} = ( XT (™, r )T (20) 

For any feasible state pair (i, j), i.e., i G [0, K] and j E [i — m, K — m], from (17) and (20), we have: 



f e -ATKr)(AT(m,r)r j<K _ m 
vl 

K-m-1 (21) 

1- Pi > z ' j = K-m 



z=i—m 



To calculate the unconditional transition probabilities, pi j, we need to average pij(m, r) over all possible 
values of m and r, i.e., 

s{i) R 

Pi,j = ^p(A m \i)^2pi,j(m,ri)(j) ri \ m (22) 

m=l 1=0 

where <p Tl \m is the probability that the rate n is chosen for a space-batch, given that it contains m packets. 
This probability is given by: 

( Pri\m / _/\^\ K^^) 

\m) 

where the sum is taken over all sets S of cardinality m, i.e., containing m distinct destination nodes. 

4) Distribution at arbitrary times, ir s : Using the PASTA property [13] of Poisson arrivals, the 
probability that at an arbitrary time in the steady-state the queue contains q = k packets is equal to 
the probability that a random arrival observes k packets in the queue. In other words, 

4 = Pr {&(*)=*} (24) 
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where ix s k is the k-th. element of 7r s , and q a (t) is the state of the queue observed by an arrival at time t. 
The right hand side of (24) can be expanded by conditioning on qd{t), the state of the queue at the most 
recent departure before t, i.e., 

k 

n s k = J2^Ua(t) = k\q d {t) = i}¥i{q d {t) = 1} (25) 

i=0 

In order to calculate Pr {qd(t) = i}, we observe that this probability can be viewed as the probability 
that the arrival at time t happens to occur during the departure state i, i.e., when the node is in state i of 
the embedded Markov chain discussed in the previous subsection. Therefore, this probability is equal to 
the expected fraction of time that the node spends in the departure state i. 

Let the random variable D{i) denote the time spent in departure state i. Of this time, T d {i) seconds 
will be spent in transmission mode, and only if i — 0, an additional / seconds will be spent in idle mode 
before entering transmission mode. Therefore: 

E[D{i)] = E[I] + E[T d {i)] = I [1 - i}+ + E[T d (i)] (26) 

where the term j- [1 — i] + is nonzero only when i = and is equal to the expected time needed for the 
queue occupancy to reach 1 packet. Then, the expected length of an interdeparture epoch is 

e \p\ = E < (\ t 1 - «r + ^woi) = + E 4e[tm (27) 

i=0 ^ ' i=0 



with 



s(i) R 

E[T d {i)} = P(A m \i)Y^ <t>n\mT(m, n) (28) 

m=l 1=0 

The probability Pr {q d (t) = i} in (25) can now be calculated as follows: 

fttoW = ,•} = «! (29) 

which can be interpreted as the fraction of total time spent in the departure state i. 

The term Pr{g a (t) = k\q d {t) = %} in (25) is the probability that an arrival during the departure state i 
observes k packets in the queue. This probability can be viewed as the fraction of arrivals in departure 
state i which observe k packets in the queue at the moment of their arrival. The expected total number 
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of arrivals in state i is given by \E[D(i)}. Of these, only one may observe k < K frames in the queue, 
provided that there are enough arrivals. Let q n+ i be the state at which the next departure will leave the 
queue. If the space-batch size for this departure is m, then the queue occupancy just before this next 
departure is q n+ i + m. In order for an arrival to have observed k < K packets in the queue, we need 
q n+ i + m > k + 1 packets. Therefore, the probability that an arrival in state i observes k packets in the 
queue, given the space-batch size m, is given by: 



R K-m 

Pr {q n +i > k + 1 - m\q n = i} = y~] <j> n \ m ^ Pi,j( m ' r i) ( 30 ) 

1=0 j=k+l-m 



AMD1 S 71- * [J2^n\m Yl PiA m ' r i) ) ) (31) 

t ^ i=0 \m=l \l=0 j=k+l-m 



The probability that an arrival in state i observes k packets in the queue is then given by: 

Pr {qM = kMt) = i} = SS.»HA. l ,)PrK»|t + l- B .| fc = i} (32) 

From (25), (29), and (32), the steady state queue occupancy distribution, iz s , for states < k < K — 1 
can be computed as shown in (31). For k = K, we have 



K-l 



7T 



K 



i=0 



Note that in (31), when k = 0, for all values of m we have 

K-m 

^ Pij(m,r) = l, 

j=k+l—m 



and (31) simplifies to: 



<=xm< <34) 



This is because in this case, during any departure state i = 0, exactly one arrival will observe packets 
in the queue with probability 1. 
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B. Performance Metrics 

Once the ir d and 7r s distributions are obtained, several performance metrics can be derived from them: 

• Blocking Probability: The probability that an arriving packet to the AP is discarded because there 
is no space for it in the transmission queue: 

• System Throughput: Number of bits that can be successfully transmitted from the AP per second: 

S = X(1- P b )L d 

• Average Queue Occupancy: The average number of packets in the queue: 



K 

£[Q] = £> 9 

q=l 



• Average Response Delay: The average delay that a packet suffers, from its entrance to the queue 
until it is transmitted, computed from the average queue occupancy by applying the Little's Law 
[13]: 

E[D] EM 



A(l - P b ) 

• Average Space-Batch Size: Average number of packets included in the transmitted space-batches: 

K I s(q) 

E[s] = XX J2 m ' Pm \i 

q=0 \m=l 

IV. Results 

In this section, the analytical model is validated through simulations, and some insights on how the 
number of antennas, number of active users, channel conditions, and traffic load impact the performance of 
a Multiuser MPT system are provided. The values of the parameters used for both the simulations and the 
analytical model are listed in Table I. A simulator of the described scenario has been built, from scratch, 
using the C++ language and based on the COST (Component Oriented Simulation Toolkit) libraries [15]. 

Figure 4 shows the behavior of different performance metrics under ideal channel conditions for two 
values of K (K = 25 and K = 100 packets) and four values of N (N = 4, N = 8, N = 16 and N = 32 
nodes) when the traffic load increases from 40 Mbps to 120 Mbps. The ideal channel conditions allow 
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Parameter 



Value 




256 bits 
64 bits 
64 bits 
8000 bits 

64 bits 
50 packets 
8 antennas 



K 
M 



{-oo, 5, 10, 15, 20, +00} 
{0, 6, 12, 18, 24} Mbits/second 



K 

T 



10 msecs 



TABLE I 

Parameters considered in Figures 4 and 5 



the AP to always transmit at the highest available transmission rate (24 Mbps), regardless of the number 
of parallel streams being transmitted. Therefore, s max can be set to its highest value (s max = M = 8) 
to maximize the number of packets that can be included in each space-batch. As a result, the system 
performance in terms of blocking probability (Figure 4(a)) and average delay (Figure 4(b)) is only affected 
by the ability of the AP to schedule large space-batches (Figure 4(c)), which in turn depends on the queue 
occupancy (Figure 4(d)) and the number of nodes sharing the aggregate traffic load. Obviously, in the 
case where there are fewer nodes than s rnax , the system performance is limited by the number of nodes. 

For K = 25, in terms of blocking probability and delay, it can be observed that as the number of nodes 
increases, the system performance improves since larger space-batches are more likely. For K = 100, 
the same behavior can be observed, although the performance gain achieved by increasing the number 
of nodes is less significant. This is because a larger queue can store more packets, and therefore, the 
probability that it contains packets directed to a higher number of destinations is also higher, which 
allows for transmission of larger space-batches even when the number of nodes is small. However, on the 
downside, a larger queue size results in a longer average delay due to increased waiting time. 

In terms of accuracy, the precision of the analytical model improves as the queue size grows. This 
means that (5) provides a better estimation of the space-batch size when the queue occupancy is higher 
and there is potentially a higher number of destinations represented in the queue. In terms of the number 
of nodes, for a fixed number of antennas (here M = 8), the model is least accurate when the number of 
nodes is exaclty equal to the number of antennas, i.e., N = M. When N < M, each transmission will 
contain as many packets as there are nodes represented in the queue. This results in a maximum loss 
of diversity in the queue immediately after a departure. However, since the largest possible space-batch 
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size is itself equal to N, for smaller N, the aforementioned maximum diversity loss is also smaller, and 
therefore, the queue is more likely to have a diverse representation of nodes after a departure. As iV 
increases, the toll of departures on the queue diversity becomes more significant until the number of 
nodes reaches N = M. For iV > M, the opposite behavior is observed, i.e., as the number of nodes 
increases, the analytical model becomes more accurate. This is because in this case the space-batch size 
is often smaller than the queue diversity, and therefore the departure has a less significant effect on the 
remaining diversity of the queue. This is why in Figure 4 the analytical curves corresponding to N = 8 
are the curves showing the highest discrepancy with their simulated counterparts. 

Figure 5 shows the behavior of the system in non-ideal and heterogeneous channel conditions. In Figures 
5(a) and 5(b), the blocking probability and the delay for N = 16 nodes and different s max values are 
plotted against different traffic loads. The 16 nodes are distributed in the coverage area of the AP, in a 
way that a first group of 5 nodes observe an average SNR equal to 25 dBs, a second group of 5 nodes 
observe an average SNR equal to 45 dBs, and finally, a third group of 6 nodes observe an average SNR 
equal to 35 dBs. As it can be seen in these figures, lower blocking probability and delay can be achieved 
by appropriately choosing the value of s max . 

In terms of blocking probability, there seems to be an optimal value of s max for every given arrival 
rate. At very low arrival rates, the packets rarely get accumulated in the queue and there is no need for 
sending large space-batches. In fact, sending large space-batches needs longer service time, which may 
increase the risk of arrivals getting blocked. This is why at low arrival rates, s max = 2 minimizes the 
blocking probability. As the arrival rate increases, the need for sending larger space-batches becomes 
more evident and s max = 4 and subsequently s max = 6 become the optimum choices. In the scenario 
considered, s max = 8 is never optimum. This is because for our particular choice of SNR values, the loss 
in terms of service time when s max = 8, which is partly caused by the lower transmission rates that have 
to be chosen for larger space-batches, surpasses the gain from the efficiency of sending more packets at 
once. 

The s max value is a trade-off between the number of packets included in each space-batch and the 
transmission rate at which the space-batch can be transmitted, as observed in Figures 5(c) and 5(d), 
respectively. When the traffic load increases, the queue occupancy grows, and the probability to schedule 
larger space-batches also increases. However, as the CSI is not used for selecting neither the number nor 
the specific destinations to which packets are sent, the transmission rate at which space-batches are sent 
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decreases as the number of spatial streams increases. Therefore, using the s max parameter, the system can 
achieve better performance by trading off the maximum number of packets transmitted with the average 
transmission rate observed. 

V. Conclusions 

We have presented a queuing model for the performance evaluation of Multiuser MPT systems. The 
analysis is built around a blind estimation of the space-batch size distribution, which is only based on the 
distribution of the traffic load between nodes, which is assumed to be homogeneous. This approximation 
allows to keep the model reasonably simple but, as the results show, also very accurate. 

In such conditions, the model is expected to be used both alone, to evaluate the impact of the number of 
nodes, number of antennas, transmission rates, etc. on the performance, or coupled with other link-layer 
mechanisms for the evaluation of more complex systems. For example, it can be easily combined with a 
model of the Distributed Coordination Function (DCF) [16] to evaluate the performance in non- saturation 
conditions of the upcoming IEEE 802.1 lac standard [17], that will support Multiuser MPT. 

In its current form, the presented model contains all the elements to characterize a Multiuser MPT 
system. However, it can be further extended in future works to cover other aspects, e.g.: 

• to consider non-uniform traffic distribution among nodes, as well as other Markov-based arrival 
processes. 

• to formulate schedulers that also consider the existing multiuser diversity (i.e., schedulers that pick 
the packets from the buffer based on the instantaneous CSI). 

• to consider transmission errors and retransmissions, in addition to multiple transmission rates. 

• to combine the model with packet fragmentation and aggregation techniques. 
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Fig. 1. Schematic model of an AP equipped with M antennas. SCH stands for scheduler. 
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Fig. 2. Frame structure of an Space-batch for a MPT system where the AP has M antennas. 




Fig. 3. Temporal evolution of the AP's queue with M = 2 antennas. The q values are the queue state (number of frames in the queue) 
just after a frame arrival or a space-batch departure. The @x refer to the packet destination. 
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Fig. 4. Performance evaluation for different N and K values in ideal channel conditions. 
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